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We report on a lattice computation of the second moment of the pion matrix element of the twist-2 
non-singlet operator corresponding to the average momentum of parton densities. We apply a fully 
non-perturbatively evaluated running renormalization constant as well as a careful extrapolation of 
our results to the continuum limit. Thus the only limitation of our final result is the quenched 
approximation. 



1 Introduction 

A reliable computation of parton distribution 
functions from first principles would be very 
important for future experiments planned 
e.g. at the LHC. The results of such cal- 
culations would offer a unique way to test 
whether QCD is indeed the correct theory for 
the strong interactions. 

In principle, the lattice regularization of- 
fers such a calculational scheme, where the 
only starting point is the QCD lagrangian, 
and indeed is able to give results for mo- 
ments of parton density distributions in the 
continuum and fully non-perturbatively, as we 
will show in this contribution. The work pre- 
sented here is summary of a series of publi- 
cations Id'H'HO'El that provide the essential in- 
gredients to reach the above ambitious aim. 

The moments of parton density distri- 
butions are related to expectation values of 
certain local operators, which are renormal- 
ized multiplicatively by applying appropriate 
renormalization factors Z(X/fi) that depend 
on the energy scale \x. This leads to con- 
sider renormalized matrix elements Og™(^) 
to be computed in a certain, at this stage not 
specified, renormalization scheme SF. If the 
energy scale [i is chosen large enough, it is to 
be expected that the scale evolution is very 
well described by perturbation theory, giving 
rise to the following definition of a renormal- 
ization group invariant matrix element: 



with g(fi) the running coupling and 
/ SF (3 2 ) = (5 2 (^)r W2h0 (2) 
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where f3(g) and 7(5) are the (3 and 
anomalous-dimension functions computed to 
a given order of perturbation theory in the 
specified scheme, i.e. here the SF scheme. 
Once we know the value of 0\^ y evaluated 
non-perturbatively, the running matrix ele- 
ment in a preferred scheme can be computed, 
for example in the MS scheme: 



0^) = 0^ v /f MS (g 2 (p)) 



(3) 



with now, of course, the /3 and 7 functions 
computed in the MS scheme. 

A non-perturbatively obtained value of 
the renormalization group invariant matrix 
element is hence of central importance. Its 
calculation has to be performed in several 
steps. The reason is that we have to cover 
a broad range of energy scales - from the 
deep perturbative to the non-perturbative re- 
gion. With the scale-dependent renormaliza- 
tion factor we can write the renor- 
malized matrix element of eq. ([j]) as 

(ttIOnsK) 



Of 



f bt (g 2 (p)), (4) 



I Tv = Or F n (M)-/ SP (5 2 (M)) 



(1) 
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with (ttIOnsI 71 ") the expectation value of the 
(non-singlet) operator under consideration in 
given states, here the pion states. 

So far, all our discussions have been in 
the continuum. However, if we think of the 
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lattice regularisation and eventual numeri- 
cal simulations to obtain non-perturbative re- 
sults, it would be convenient to compute the 
renormalized matrix element at only one con- 
venient (i.e. small hadronic) scale fio- We 
therefore rewrite the r.h.s. of eq. (0) as: 



/ SF Gf(M)) (5) 



=<t(/x/mo,9(a0) 



where we introduce the step scaling function 
cr(/i//i , g{^)), which describes the evolution 
of the renormalization factor from a scale [io 
to a scale /i. The advantage of concentrating 
on the step scaling function instead of the 
renormalization factor itself, is that the step 
scaling function is well defined in the contin- 
uum and hence suitable for eventual contin- 
uum extrapolations of lattice results. 

We finally write the r.h.s. of eq. (Q) as 

Osf(w>) aQi/toM) ■ / SF (ff 2 (M)) (6) 

= S« v v (po) 

with Og™(/io) the renormalized matrix ele- 
ment, which is to be computed only once 
at a scale /iq and the (ultraviolet) invariant 
step scaling function 6}jJy(/io), which still 
depends on the infrared scale /io- The fol- 
lowing sections are devoted to a description 
of how these two basic ingredients can be 
reliably computed on the lattice using non- 
perturbative methods, i.e. numerical simula- 
tions. 

2 The renormalization group 
invariant step scaling function 

Let us start this section by disclosing what is 
hidden behind the fictitious SF scheme men- 
tioned in the introduction. SF stands for 
Schrodinger functional and denotes a finite 
physical volume, V — L 3 ■ T, renormaliza- 
tion scheme where the energy scale fi is iden- 
tified with the inverse spatial length of the 
box itself, e.g. /i = l/L. The peculiarity of 
the Schrodinger functional set-up is that fixed 



boundary conditions in time xq are imposed 
with classical fields at the time boundaries at 
time xo = and xq = T. For a more detailed 
discussion we refer the reader to 0. 

To discuss the renormalization of oper- 
ators related to moments of parton distri- 
bution functions, we first have to provide a 
renormalization condition. Denoting by \SF) 
a classical SF state, i.e. a classical quark field 
at a time boundary with external momentum 
p, the renormalization condition that we will 
use reads 



(SF\O r 



\SF) = (SF\O tICC \SF) 



(7) 

The relation between the expectation value of 
the bare operator and the renormalized one is 
established through a scale-dependent renor- 
malization constant: 



R {n) = Z-\l/fx)O haIe (l/L) . (8) 

In perturbation theory, on the 1-loop level, 
we have Z(l//i) = 1 - f^i) [ 7 (°)ln(/i) + B ] 
with the anomalous dimension and B$ 
the constant part. Up to this point, the dis- 
cussion is given solely in the continuum where 
the SF renormalization scheme is a perfectly 
acceptable one. Different schemes such as the 
MS scheme can be related to the SF scheme 
as usual in perturbation theory. 

If we are interested, however, in a non- 
perturbative calculation, we have to detour 
for a short time (which means, however, a 
substantial computer time) to a finite lattice 
with non-zero lattice spacing a that allows 
for numerical simulations. A lattice repre- 
sentation of the twist-2 non-singlet operator, 
which is the only case we are considering here, 
is given by 

1 - - r 3 

12 (x) = -i>(x) 1{1 D 2} —1>(x) , (9) 

where is the covariant derivative and the 
bracket around indices means symmetriza- 
tion. The operator is probed by boundary 
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W-) = £ <e ip(y - z) Oi 2 (x)C(y) 72 r 3 C(z)) 

n ^ — ' 



quark fields ( and which reside at a:o = 
and a correlation function is constructed 

a 

x,y,z 

(10) 

with p the spatial 3-momentum. 

To take into account the effects of the 
boundary fields, we also consider the bound- 
ary operators defined at the time boundaries 
xq = and xq = T: 



from which we construct the correlation func- 
tion: 

fi = -(O O T ). (12) 

The boundary wave-function contribution 
can then be taken out by considering the ra- 
tio fo ia (vo/a)/V7i- 

There are several physical scales in our 
problem which all have to be given in units 
of L, which is the only scale that is to be 
changed. Therefore, in order to arrive at a 
definition of the renormalization constant, we 
have to specify the spatial momentum p, the 
quark mass m q and the time xq where we 
read off the expectation value of the operator 
between SF states from the correlation func- 
tion. The final physical result does, of course, 
not depend on our choice of these quantities, 
but we choose them solely for convenience. 
In particular we select 

m q = 0, x =T/A, (13) 

P = (P1,P2,P3) = (Pi = 2vr/L,0,0) . 

The choice of a zero quark mass results in 
using a massless renormalization scheme and 
the choice of the smallest available momen- 
tum on the lattice minimizes lattice artefacts. 
With the above choice, it is indeed only the 
physical box length L (assuming T = L) that 
we identify with the inverse scale, which is 
varied in the problem. 



We are now in a position to give the 
precise definition of the renormalization con- 
stant 



Z(L) = Z{L)/Z 1 (L) , (14) 



with 



Z(L) = f 0l2 (T/4)/fS™(T/4), 

Z^L) = yfimiylft~(L) , (15) 

where we divide by the corresponding tree 
level expression as required by the renormal- 
ization condition, eq. (Q). Instead of comput- 
ing the Z-factors, we concentrate on the step 
scaling functions 



Z(2L) 



°7i = 



Z X {2L) 
Zi(L) ' 



Z(2L) 



Z(L) ' 
(16) 

because, in contrast to the Z-factors, the step 
scaling functions have a well-defined contin- 
uum limit. The strategy is now to compute 
the step scaling functions at various values 
of the lattice spacing while keeping fixed the 
conditions in eqs. (|l^) and the physical scale 
ji = L^ 1 (determined by the running cou- 
pling g(fj,)) and to extrapolate the results 
thus obtained to a = 0. 

It is one of the basic ingredients and char- 
acteristics of our work that almost all simula- 
tion results at non-zero lattice spacings have 
been obtained by employing the standard 
Wilson action and the non-perturbatively im- 
proved clover action. Since these two formu- 
lations lead to different lattice artefacts, it is 
a very crucial test of our results that their 
continuum extrapolations give consistent re- 
sults. That this is indeed the case is demon- 
strated in fig. |l|. It shows that for the two step 
scaling functions Uf t and a 2 the continuum 
limit of both discretizations agree within the 
error bars. We note that in the case of a 2 a 
quadratic extrapolation in the lattice spacing 
a is necessary while for af t a linear extrap- 
olation is sufficient. After checking that a 
similar behaviour is found at all values of the 
coupling we have simulated, we performed 
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Figure 1. Continuum extrapolation of the step 
scaling functions erg and cry. performed separately 
for the Wilson action (circles and dotted lines) and 
the non-perturbatively improved action (squares and 
full lines) at a fixed value of the running coupling 
g 2 = 2.45. 
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Figure 2. Constraint continuum extrapolation of a % 
(notation as in fig. 



constraint fits, demanding that the contin- 
uum value of the step scaling functions be 
the same from both actions. A summary of 
our results for a 2 is shown in fig. |[ 

At this point, having performed all nec- 
essary continuum extrapolations, we end our 
detour on a lattice with non-zero lattice spac- 
ing and come back to the discussion in the 
continuum. With the results on the step scal- 
ing function, extrapolated to the continuum 
limit, which were obtained at 9 values of the 
running coupling constant, we can now com- 
pute the (ultraviolet) invariant step scaling 
function 

®inv(Mo) = er(/i/Mo,.9 2 (W))) ■ fifiv)) , 

(17) 

which still depends on the infrared scale /iq. 
This scale dependence will only be cancelled 
when multiplying with the matrix element, 
renormalized at the scale (j,q. The function 
f{g 2 (^)) is the same as in eq. (||) and the (3 



and 7 functions are taken up to 3 loops in the 
SF schemeR 

In fig. H we show 6{^(^ ) as a func- 
tion of fi/ fJ.Q- For large enough energy scales, 
n/Ho > 100, 6{^,(^o) does not change 
within the errors and we can determine a 
value for it by fitting the last, say, 4 points to 
a constant. Although there still is a scheme 
dependence in the invariant step scaling func- 
tion through the remaining dependence on 
fio, the value should be independent of the 
choice of coupling used in the analysis. This 
is nicely illustrated in fig. |], where the choices 
of g(L/A) and g(L) give consistent values for 
the invariant step scaling function. 

We can therefore now give the first piece 
of information for the invariant matrix ele- 



a For the 7 function we have taken an effective 3-loop 
parametrization as obtained by fitting our data to an 
effective 3-loop form a. 
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Figure 3. The values of Sj^(/io) f° r two choices of 
the running coupling. 



ment itself, as needed in eq. (Q), and quote 



Si U N V vM = 1.11(4) • 



(18) 



3 Renormalized matrix element 



As a next step we have to compute the renor- 
malized matrix element itself. Again, we will 
always use the SF scheme and we remark that 
this is the first time that a calculation of a 
2-quark matrix element is attempted in this 
set-up. We first tried to compute the matrix 
element (within pion states) of the operator 
of eq. (S). However, since this operator needs 
a non-vanishing momentum, we found, simi- 
lar as in that the signal is very noisy. We 
therefore decided to switch to the operator 



' « i 3 
7o Do - g ^2lk D k 

k=l 



O 00 (x) = ip(x) 70 Do — / J7fe Dk -g-^(x) 

(19) 

which has the advantage that it can be com- 
puted at zero momentum. Taking the bound- 
ary operators of eq. (|ll|) we construct a cor- 
relation function 



(20) 



which again is to be normalized by fx, 
eq. (|l2|), to take out the boundary wave- 
function contributions. Performing a transfer 
matrix decomposition, we find that for large 
enough values of xq and staying far enough 
from both boundaries 



h 

/m(zo) 



P 2 ^ 



MO, 



oo W 



(21) 



Assuming that there is a plateau region where 
/m(x )//i = const = (7r|0 ok), and in 
which the first excited state gives essentially 
no contributions, we obtain the physical ma- 
trix element (x) after a suitable normaliza- 
tion (see □) : 



— MOook) 



(22) 



In order to extract the matrix element of 
eq. ( |22l) we have chosen lattices with T = 
3 fm and followed the correlation function 
/m(#o) up to a distance of 1 fm in time di- 
rection. At this distance we are sure that 
we project on the pion states as an inspec- 
tion of the pseudoscalar and axial-vector cor- 
relation functions (from which we also ex- 
tracted the pion masses) showed. Indeed, for 
1 fm < xq < 2 fm the correlation function 
exhibits a plateau behaviour as can be seen 
in fig. 2 of ref. la. 

Once we have the bare matrix element 
we need to renormalize it. To this end we 
computed Z(1//j,q) with fj,^ 1 = 1.436ro, ps 
0.5 fm. We repeated such a calculation for 
various lattice sizes, choosing the values of 
such that L max = Pq 1 is kept fixed. In- 
terpolating the numerical simulation data we 
obtain in this way Z(l/ '/j,q) in a range of lat- 
tice spacing 0.05 < a < 0.1, i.e. the range 
of a where the bare matrix element itself has 
been computed. 

This allowed us to apply the renormal- 
ization factor at exactly the same values of 
fj where the matrix element has been com- 
puted. We show the matrix element renor- 
malized with £>i2(l/Mo) in fig- f§ Again all 
calculations have been performed with two 
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Figure 4. Constrained continuum extrapolation of 
the renormalized matrix element 



different discretizations, and we find quite 
convincingly the same continuum limit. 

It might, however, not have escaped the 
reader's attention that we have used the 
wrong renormalization factor, namely Z\2, 
for rcnormalizing the operator Oqq{x) of 
cq. (Jig)) - The continuum extrapolation of the 
such renormalized operator needs a correc- 
tion factor. In ref. □ we have demonstrated 
that this correction factor can be taken from 
perturbation theory and amounts to a shift of 
the continuum renormalized matrix element 
by a few per cent. Taking the correction fac- 
tor to be the same as found in i, we can finally 
give our main result: 



O 



INV 



0.222(24) 



(23) 



4 Conclusion 



In this contribution we have demonstrated 
how we can compute on the lattice, in a fully 
non-perturbative fashion, a renormalization 
group invariant matrix element OJ^v f° r ^ ne 
second moment of parton distributions of the 
twist-2 non-singlet operator in the pion. A 
preliminary value for Of^v * s gi ven m eq. ( p3| ) 
and can be used now as integration constant 



to obtain the renormalized matrix element 
at any scale in the preferred renormalization 
scheme. 

We want to emphasize that the value of 
the renormalized matrix element is shifted 
substantially from a value of about (x)(a = 
0.093) = 0.30 at /3 = 6 to (x)(a = 0) = 0.2 
in the continuum limit. Hence we experience 
strong lattice artefacts in the renormalized 
matrix element. Still, when the matrix ele- 
ment is run in the MS scheme to a scale of 
fx = 2.4 GeV, we find (x)(p = 2.4GeV) w 0.3. 
Thus we find that, by a conspiration of two 
effects, our result agrees with the number 
quoted in the pioneering work of There- 
fore, the fact that within the quenched ap- 
proximation, used here exclusively, the num- 
ber from the lattice simulations is higher than 
the experimental value persists. 
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